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Abstract: Anisakis spp. (Nematoda: Anisakidae) parasitize a wide range of marine animals, mammals serving as the de- 
finitive host and different fish species as intermediate or paratenic hosts. In this study, 1 8 fish species were investigated for 
Anisakis infection. Katsuwonus pelamis, Euthynnus affinis, Caranx sp., and Auxis thazard were infected with high preva- 
lence of Anisakis type I, while Cephalopholis cyanostigma and Rastrelliger kanagurta revealed low prevalence. The mean 
intensity of Anisakis larvae in K. pelamis and A. thazard was 49.7 and 5.6, respectively. A total of 73 Anisakis type I larvae 
collected from K. pelamis and A. thazard were all identified as Anisakis typica by PCR-RFLP analysis. Five specimens of 
Anisakis from K. pelamis and 1 5 specimens from A. thazard were sequenced using ITS1 -5.8S-ITS2 region and 6 speci- 
mens from A. thazard and 4 specimens from K. pelamis were sequenced in mtDNA cox2 region. Alignments of the sam- 
ples in the ITS region showed 2 patterns of nucleotides. The first pattern (genotype) of Anisakis from A. thazard had 100% 
similarity with adult A. typica from dolphins from USA, whereas the second genotype from A. thazard and K. pelamis had 
4 base pairs different in ITS1 region with adult A. typica from USA. In the mtDNA cox2 regions, Anisakis type I specimens 
from A. thazard and K. pelamis showed similarity range from 94% to 99% with A. typica AB517571/DQ1 16427. The dif- 
ference of 4 bp nucleotides in ITS1 regions and divergence into 2 subgroups in mtDNA cox2 indicating the existence of A. 
typica sibling species in the Makassar Strait. 
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INTRODUCTION 

The Indonesian archipelago, which comprises of more than 
17 thousands islands and coastal lines length of about 104,000 
km, provides great natural resources for Indonesia economic 
growth in many regions, particularly the fishery sector. Total 
production from fisheries continues to grow from 8.2 million 
ton in 2007 to 12.4 million ton in 2011 [1]. Indonesia is also 
recognized as a center of marine organism diversity and is in- 
habited with abundant fish species. Fishery products, which 
are either traded as export commodities or for local consump- 
tion, have to be free from any zoonotic parasites, such as ani- 
sakid nematodes. Food safety has recently become a great con- 
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cern for consumers. Zoonotic parasites in fish products are 
mainly caused by helminths, which utilize fish as the interme- 
diate host, and marine mammals such as dolphins or whales 
as the final host [2-4]. Some helminths are transmittable and 
able to survive in the human digestive tract after consuming 
raw fish infected by larval helminths, causing significant clini- 
cal diseases, as well as allergic reactions [5-7]. Symptoms of 
anisakiasis include epigastric pain, nausea, vomiting, and diar- 
rhea [8]. Among the helminth parasites, larval nematodes of 
the genus Anisakis are commonly found in the musculature or 
digestive tracts of many species of marine fish. Anisakis simplex 
is the most well-known zoonotic nematode which has been 
reported to cause anisakiasis of humans in many countries in 
Europe and Asia [8]. Though reported cases of anisakiasis in 
Indonesia are very rare, a study conducted by Uga et al. [9] us- 
ing a seroepidemiological approach of inhabitants in East Java 
revealed that about 11% of samples were positive for Anisakis 
infection. 

In previous studies, Anisakis spp. could only be categorized 
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morphologically into Anisakis type I and type II, in which the 
former has a longer ventriculus and a mucron, while the latter 
has short ventriculus and no mucron [10]. Identification to the 
species level by microscopic examinations is usually unreliable 
due to undeveloped morphological characteristics of larval stage 
nematodes. However, accurate identification of Anisakis nema- 
todes larvae is required for precise diagnosis of Anisakis infec- 
tions in humans and fish and to improve food safety. In addi- 
tion, they can be used as biological indicators in the study of 
stock discrimination of migratory fish [11-15]. Recent studies 
showed that molecular diagnostic techniques could be used to 
identify Anisakis to a species level. PCR-RFLP has been widely 
used to identify Anisakis spp. in different fish species [16-20], 
Molecularly, the previous Anisakis type I is known to consist of 
6 species (A. ziphidarum, A, nascettii, A. typica, and 3 sibling spe- 
cies of A. simplex complex, namely A. simplex (sensu stricto) 
(s.s.), A. pegreffii, and A. simplex C, whereas type II consists of 3 
species (A. paggiae, A. brevispiculata, and A. physeteris) [21,22], 
Recent studies in Japan showed that L3, L4, and adult of A. 
simplex (s.s) and A. pegreffii could be distinguished morpho- 
logically based on the ventricululs length, in which the former 
has a longer (0.90-1.50 mm) ventriculus than the latter (0.50- 
0.78 mm) [23]. 

Parasitological research on Anisakis spp. in Indonesia is rela- 
tively scarce. Research on Anisakis was previously reported from 
Seribu Islands, Jakarta from 3 species of fish: Rastrelliger kana- 
gurta, Decapteivs russelii, and Sardinella sirm [24]. The study has 
found 2 types of Anisakis, i.e., Anisakis type I and Terranova type 
B, where Anisakis type I predominated. A recent study conduct- 
ed on several fish species of Balinese and Javanese waters has 
found 3 species of Anisakis; A. typica, Anisakis sp. 1, and Anisa- 
kis sp. 2 [25]. Another study of marine fish in Indonesia in 
Southern Coast of Kulon Progo, Yogyakarta, has also found 5 
out of 11 fish species examined harbored Anisakis spp. [26]. 
However, from most studies conducted in Indonesia, except 
for that conducted by Palm et al. [25], identification of the An- 
isakis larva was solely based on morphology, making it unreli- 
able for species identification. 

The aim of the present study was to investigate the occur- 
rence of Anisakis infection from some marine fish in the South- 
em Makassar Strait, and characterize them to species level us- 
ing PCR-RFLP genetic analysis, the molecular keys described 
by D'Amelio et al. [27] and Pontes et al. [17], and sequencing 
of ITS-5.8S and mitochondrial cytochrome c oxidase subunit 
II (mtDNA cox!) regions. 



MATERIALS AND METHODS 

Parasitological examinations 

A total of 220 fish representing 18 species and 10 families 
were investigated for Anisakis larvae infection. Fish species, to- 
tal length, and number of fish examined are given in Table 1. 
Fish were purchased from local markets located in Makassar, 
Takalar, and Barm regencies (Fig. 1). The fish were transported 
to the Laboratory of Fish Parasites and Diseases, Hasanuddin 
University for parasitological examinations. The length of fish 
was measured, body cavity was opened, and internal organs 
were placed on Petri dishes and examined for the presence of 
Anisakis larvae. Anisakis larvae were distinguished from other 
anisakids following of the identification protocols of Ander- 
son [28]. Anisakis spp. were then categorized into type I or type 
II. Anisakis type I was characterized by the presence of boring 
tooth at the anterior end, ventriculus, and mucron at the pos- 
terior end, whereas type II has a boring tooth, ventriculus, and 
no mucron at the posterior end. Anisakid larvae were isolated 
from the visceral surface and body cavity of the fish. The larvae 
were observed under light and dissection microscope for mor- 
phological identification. All Anisakis larvae found were count- 
ed. Since the number of Anisakis larvae found from fish other 
than skipjack and frigate tuna was very few (mostly uninfect- 
ed), Anisakis spp. from these fish were only morphologically 
identified into Anisakis type I or type II. Only larval Anisahs spp. 
isolated from skipjack tuna and frigate tuna were used for fur- 
ther molecular characterization. 

The prevalence is defined as the percentage of fish infected 
by Anisakis larvae. The mean intensity is defined as the total 
number of Anisakis species found divided by the total number 
offish infected [29]. 

DNA isolation and amplification 

Anisahs spp. were fixed and stored in 70% ethanol. Genom- 
ic DNA from individual worms was extracted using a QIAamp 
DNA Mini Kit (Qiagen Inc., Hilden, Germany) following the 
manufacturer's tissue protocol. DNA was eluted with milliQ 
and stored at -20°C before subsequent PCR amplification. In- 
ternal transcribed spacer (ITS) and mitochondrial cytochrome 
c oxidase subunit II (mtDNA cox2) regions were used for am- 
plification and sequencing. The entire rDNA regions compris- 
ing ITS1, 5.8S, and LTS2 were amplified using previously de- 
scribed primers NC5 (5 -GTAGGTGAACCTGCGGAAGGAT- 
CATT-3') and NC2 (5'-TTAGTTTCTTTTCCTCCGCT-3) [30], 



Anshary et al.: Molecular identification oMn/sa/c/s from marine fish, Indonesia 11 



Table 1 . Fish species, number examined, prevalence, and intensity of Anisakis type I infection 









Number of fish 


(%) 


Intensity 


Fish species 


Locality/Size 


Date 


Examined 


Infected 


Mean intensity 
(range) 


Caesionidae 














Fusilier Caesio sp. 


Barru/25-32 cm 


79/10 


9 


0 


0 




Carangidae 














Giant trevally Caranx sp. 


Makassar/64-67.5 cm 


26/8/10 


4 


3 


75 


1 (1) 


Indian scad Decapterus russelii 


Makassar/23 -36.5 cm 


26/8/10 


20 


0 


0 




Longnose trevally Carangoides sp. 


Barru/25-30 cm 


79/10 


g 


0 


0 




Clupeidae 














Goldstripe sardinella Sardinella sp. 


Makassar 


16/9/10 


41 


o 


0 




Lutjanidae 














Snapper Lutjanus sp. 


Barru/22-29 cm 


79/10 


g 


o 


0 




Priacanthidae 














Priacanthus sp. 


Barru/28 cm 


79/10 


i 

1 


U 


0 




Scaridae 














Parrot fish Scarus sp. 


Barru/28 cm 


79/10 


4 
1 


V 


0 




Scornbridae 














Frigate tuna Auxis thazard 


Takalar/19-25 cm 


26/8/10 


1 2 


0 


0 






Makassar/33-41 cm 


16,17,22/9/10 


on 
oU 


1 4 


46.7 


5.6 (1-19) 


Indian mackerel Rastrelliger kanagurta 


Makassar/20-30.5 cm 


26/8/10 


1 9 
I <i 


n 
u 


0 






Takalar/20.2-27 cm 


22/9/10 


on 


I 


5 


1 


Mackerel tuna Euthynnus affinis 


Makassar/52-56 cm 


9/8/10 


Q 

o 


o 


66.7 


1 (1) 




Barru/51-52 cm 


79/10 


o 
o 


n 
u 


0 




Narrow-barred Spanish Scomberomorus 


Makassar/38.5-45 cm 


20/9/10 


4 


0 


0 




commersoni, 














Skipjack tuna Katsuwonus pelamis 


Makassar/35-60 cm 


26/8/10 


13 


12 


92.3 


49.7(1-175) 


Serranidae 














Grouper Cephalopholis cyanostigma 


Makassar/22-25 cm 


26/8/10 


8 


1 


12.5 


1 




Makassar/20-24 cm 


16/9/10 


10 


0 


0 




Grouper Epinephelus fuscoguttatus 


Barru/26-29 cm 


79/10 


2 


0 


0 




Siganidae 














Barhead spinefoot Siganus virgatus 


Barru 


79/10 


1 


0 


0 




Goldlined spinefoot Siganus guttatus 


Makassar/23-32 cm 


26/8/10 


5 


0 


0 




Sphyraenidae 














Barracuda Sphyraena sp. 


Barru/37-49 cm 


79/10 


4 


0 


0 





All PCR was performed in 20 ul which contain approximately 
dNTP 0.2 mM, primers 0.8 uM, Taq polymerase 0.02 U/pl and 
10 x buffer PCR IX, and 2 pi samples. Milli-Q was added to 
achieve the total PCR volume. Each PCR reaction was per- 
formed in a thermocycler iCycler (Bio-Rad, Hercules, Califor- 
nia, USA) under the following conditions: after initial denatur- 
ation at 95°C for 15 min, 30 cycles of 94°C for 1 min (denatur- 
ation), 55°C for 1 min (annealing), 72°C for 1 min (extension), 
followed by a final extension at 72°C for 5 min. The mtDNA 
cox2 gene was amplified using the primers 210 (5'-CAC- 
CMCTCTTAAMTTATC-3') and 211 (5'-TITCTAGTrATATAG- 
ATTGRT-TYAT-3') [31]. The PCR mixture was denatured at 94°C 
for 3 min, followed by 35 cycles at 94°C for 30 sec, 46°C for 1 
min, 72°C for 1 min and 30 sec, followed by post-amplification 



at 72°C for 10 min. The PCR products obtained were visualized 
in an SYBR green stained 1.5% agarose gel. 

PCR-RFLP 

PCR products of about 965 bp amplified with primers NC5 
and NC2 were used for PCR-RFLP analysis to identify Anisakis 
spp. following D'Amelio et al. [27] and Pontes et al. [17]. Three 
individual restriction enzymes [Taq I, Hinf I, and C/o 1) were 
used. The PCR products were digested following the manufac- 
turer's recommendation. Briefly, amplicons of 8 pi were mixed 
with 10 x reaction buffer, 0.5% BSA (only for Taq I), and di- 
gested with resttiction enzymes Taq I (10 U/pl, Takara) at 65°C 
for 3-4 hr, and with Hinfl (10 U/pl, Roche) and C/o I (10 U/ul, 
Roche) at 37°C for 3-4 hr. Milli Q was added to reach a final 
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Fig. 1 . Sampling sites of marine 
fish at the Southern Makassar 
Strait: Barru, Makassar, and 
Takalar Regencies, Indonesia. 



volume of 20 pi. The digested samples were then separated by 
electrophoresis using 1.5% agarose gel at 100 V for 40 min, 
stained with SYBR green, and photographed. Their size was es- 
timated using 100 bp ladder marker (Takara). 

Sequencing 

Twenty specimens, 5 from K. pekmis and 15 from A. thazard 
were molecularly identified using PCR-sequencing in ITS1- 
5.8S-ITS2. In addition, 10 Anisakis samples (6 from A. thazard 
and 4 from K. pekmis) were sequenced in mtDNA cox2 region. 
PCR products were purified using a PCR purification kit (Qia- 
gen) and used directly in sequencing reactions. A 100 bp or 1 
kb ladder marker (Takara) was used to estimate the size of PCR 
products. Afterwards, a total volume of 14 pi containing 6.4 
pmol primer and 10 to 40 ng DNA was prepared and sent to 
Operon Biotechnologies Company (Tokyo, Japan) for sequenc- 
ing. Milli-Q was added when necessary for DNA dilution to 
meet the concentration of DNA required. Both spacers (ITS1 
and ITS2) and the 5.8S gene were sequenced in both directions 
from each PCR product, using the same primers as above (NC5 
and NC2), NC13 (forward; 5'-ATCGATGAAGAACGCAGC-3'), 
NC13R (reverse; 5 '-GCTGCGTrCTTCATCGAT-3 '), and XZ1R 
(reverse; 5'-GGAATGAACCCGAIGGCGCAAT-3'). Both forward 
(primer 210) and reverse (primer 211) directions of mtDNA 
cox2 region was sequenced using the same primer as used for 
PCR amplification. 



Alignment and phylogenetic analysis 

The forward and reverse sequences of ITS (ITS1, 5.8S, and 
ITS2) and mtDNA cox2 regions were assembled and edited us- 
ing Bioedit Allignment Sequence Editor Ver. 7.0.5.3. They were 
compared manually with the original chromatograms when 
necessary. The obtained sequences were aligned with previous- 
ly characterized sequences of Anisakis spp. registered in Gen- 
Bank, using CLUSTAL X Version 2.1 Multiple Sequence Align- 
ments [32]. Phylogenetic and molecular evolutionary analyses 
were conducted using MEGA version 5 [33]. Maximum likeli- 
hood tree was constructed for ITS-5.8S using Pseudoterranova 
decipiens as an outgroup and neighbour-joining tree for mtD- 
NA cox2 using H. reliquens as an outgroup. ITS-5.8S and mtD- 
NA cox2 gene sequences were deposited in GenBank under ac- 
cession no. KC928261 to KC928272. 

RESULTS 

Prevalence and intensity of Anisakis larvae 

Morphologically, all Anisakis larvae found were identified as 
Anisakis type I. Among the fish examined, K. pekmis, E. affinis, 
Caranx sp., and A. thazard were infected at a high prevalence 
and high intensity, while other fish (C. cyanostigma and R. ka- 
nagurta) were infected with a low prevalence and low intensity 
(Table 1). The remainders were not infected. The parasites were 
mainly found on the surface of the viscera such as the liver and 
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Table 2. PCR-RFLP patterns Anisakis species [36] and the present samples using Taq I, Hinf I or Cfo l/Hha I restriction endonucleases 



Anisakis species 




Enzymes 




Taql 


Hinf I 


Cfo l/Hhal 


Anisakis simplex s.s. 


430, 400, 100 


620, 250, 80 


550, 430 


A. pegreffii 


400, 320, 150 


370, 300, 250 


550, 430 


A. ziphidarum 


330, 300, 140 


370, 320, 290 


550, 430 


A. typica 


400, 350 


620, 350 


320, 240, 180, 160 


Present sample 


400, 350 


620, 350 


320, 240, 180, 160 



intestines, and no Anisakis larvae were found in the muscle. 
For identification, the parasites were cleaned and the sheath 
removed. Anisakis larvae can be distinguished from other ani- 
sakid larvae such as members of Pseudoterranova, Hysterothylaci- 
um, and Contracaecum based on the shape of the ventriculus, 
which is clearly visible under a stereomicroscope. The preva- 
lence of Anisakis type I in K. pelamis was 92.3%, and mean in- 
tensity was 49.7 parasites/fish. The highest number was 175 
parasites in an individual fish. The prevalence of Anisakis type I 
larvae in A. thazard (33-41 cm in total length) reached 46.7% 
with a mean intensity of infection of 5.6 parasites/fish. Howev- 
er, none of A. thazard of smaller size (19-25 cm in total length) 
were infected out of 12 fish examined. The other fish with a 
high prevalence of Anisakis type I were Caranx sp. (75%) and E. 
affinis (66.7%) (Table 1). 

PCR-RFLP patterns 

Amplification of entire ITS and 5.8S regions of all specimens 
of Anisakis produced a PCR product of about 960 bp. The PCR 
products were digested using 3 different restriction enzymes, 
Taq I, Hinf I, and Cfo I. In PCR-RFLP, all specimens digested with 
Taq I, Hinf I, and Cfo I indicated that the samples belong to A. 
typica (Table 2). Based on the RFLP analyses of 73 Anisakis type 
I specimens from K pelamis (40 specimens) and A. thazard (33 
specimens) were all (100%) identified as A. typica (Table 2). 



pairs in ITS1 region. The first nucleotide pattern was only re- 
corded from A. thazard, whereas the second one was found 
from both K. pelamis and A. thazard. The first pattern (geno- 
type) showed 100% similarity with adult A. typica reported 
from dolphins in USA and high similarity with the Brazilian A. 
typica from dolphins, which only differed in 3 deletions in 
rrsi. Whereas the second pattern (genotype) has 4 base pairs 
difference in ITS 1 with A. typica from USA, but 100% similarity 
with Indonesian A. typica EU346093 from fish Aims rochei ro- 
chei, 4 base pairs different with A. typica EU346092 and 2 base 
pairs different with A. typica EU346091. A phylogenetic tree 
using maximum likelihood showed that A. typica found in the 
present study were in the same clade with other A. typica pub- 
lished in GenBank (Figs. 2, 3). Sequences of 10 samples of An- 
isakis using the primer in mtDNA coxl region produced about 
600 bp nucleotides. Pair distances of the alignment of mtDNA 
coxl showed 94-99% similarity of present samples with A. typ- 
ica AB517571 from Scomber japonicus, 93-98% with A. typica 
AB517572 from S. japonicus, and 94-100% with adult A. typica 
DQ116427 from dolphins (Table 3). The phylogenetic tree of 
mtDNA coxl region showed that all samples were in the same 
duster with A. typica but produced broad divergence consisting 
of 2 subgroups (Fig. 4). The first subgroup showed 96% to 
100% similarity with the known A. typica, whereas the second 
one has 93-95% similarity with the nematode (Table 3). 



Sequencing of entire ITS region and mtDNA cox2 

Sequencing of entire ITS and 5.8 S regions was performed 
for 5 samples of Anisakis type I from K. pelamis and 15 from A. 
thazard. Using the primer in the ITS and 5.8S region approxi- 
mately 950 bp nucleotides were generated. Nucleotide sequenc- 
es from all specimens were analyzed using the software Bioed- 
it, and were manually compared with chromatogram when 
necessary. No variation in the nucleotide sequences was found 
in Anisakis from K. pelamis, whereas samples from A. thazard 
showed 2 nucleotide sequence patterns. They differed in 4 base 



DISCUSSION 

The present study provides molecular identification of Ani- 
sakis from K. pelamis and A. thazard using PCR-RFLP and se- 
quencing of ITS-5.8S and mtDNA coxl regions. This is the first 
record of molecular identification of Anisakis type I from fish 
of eastern part of Indonesia. The first molecular identification 
of Anisakis, which consisted of 3 different genotypes, namely, 
A. typica, Anisakis sp. 1, and Anisakis sp. 2, was reported from 
Balinese and Javanese waters [25]. In the present study, based 
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AY826724 
AB479120 
KC9282 61 
KC 92 82 62 
EU346093 
EU34 6092 
EU346091 



AY826724 
AB479120 
KC928261 
KC928262 
EU34 6093 
EU346092 
EU346091 



AY826724 
AB479120 
KC928261 
KC9282S2 
EU346093 
EU346092 
EU34 6091 



AY826724 
AB479120 
KC 928261 
KC928262 
EU346093 
EU346092 
EU346091 



AY826724 



10 

CAACGTGCAT ACCGCCCATT 
CAACGTSCAT ACCGCCCATT 
CAACGTGCAT ACCGCCCATT 
CAACGTGCAT ACCGCCCATT 
CAACGTGCAT ACCGCCCATT 
CAACGTGCAT ACCGCCCATT 
CAACGTGCAT ACCGCCCATT 



2 0 



30 



4 J 



50 



'. 0 CATGTTGT 

r . C CATGTTGT 

r . 0 CATGTTGT 

r . a CATGTTGT 

'. i CATGTTGT 

'. G CATGTTGT 

'. G CATGTTGT 



TGTGAGCCGC 
TGTGAGCCGC 
TGTGAGCCGC 
TGTGAGCCGC 
TGTGAGCCGC 
TGTGAGCCGC 
TGTGAGCCGC 



ACGGAAACTC 
ACGGAAACTC 
ACGGAAACTC 
ACGGAAACTC 
ACGGAAACTC 
ACGGAAACTC 
ACGGAAACTC 



GTACACGTTT 
GTACACGTTT 
GTACACGTTT 
GTACACGTTT 
GTACACGTTT 
GTACACGTTT 
GTACACGTTT 



■ 1 . 



■ I 



- I 



. I. 



■ I 



.1 ■ 



■ I 



■ I 



■ I 



70 80 90 100 110 120 

GTGGTGGTGA TAGCCGTCTG CTGTGCGTTC GTTGGGCAGA CAATGGCTTA CGAGTGGcBt 

GTGGTGGTGA TAGCCGTCTG CTGTGCGTTC GTTGGGCAGA CAATGGCTTA CGAGTGGC - T 

GTGGTGGTGA TAGCCGTCTG CTGTGCGTTC GTTGGGCAGA CAATGGCTTA CGAGTGGC - T 

GTGGTGGTGA TAGCCGTCTG CTGTGCGTTC GTTGGGCAGA CAATGGCTTA CGAGTGGC - T 

GTGGTGGTGA TAGCCGTCTG CTGTGCGTTC GTTGGGCAGA CAATGGCTTA CGAGTGGC - T 

GTGGTGGTGA TAGCCGTCTG CTGTGCGTTC GTTGGGCAGA CAATGGCTTA CGAGTGGC - 1 

GTGGTGGTGA TAGCCGTCTG CTGTGCGTTC GTTGGGCAGA CAATGGCTTA CGAGTGGC - 1 



I | | | 1 | | | | | . 

130 140 150 160 170 

GTGC C SCTTG TTGAACAACG GTGACCAATT TGGC T 3TCTA CGCCGTATCT AGK 

GTGC - SCTTG TTGAACAACG GTGACCAATT TGGC - 3TCTA CGCCGTATCT AGK 

GTGC - GCTTG TTGAACAACG GTGACCAATT TGGC - 3TCTA CGCCGTATCT AGi 

GTGC - SCTTG TTGAACAACG GTGACCAATT TGGC - 3TCTA CGCCGTATCT AGK 

GTGC - GCTTG TTGAACAACG GTGACCAATT TGGC - 3TCTA CGCCGTATCT AG< 

GTGC - GCTTG TTGAACAACG GTGACCAATT TGGC - 3TCTA CGCCGTATCT AGI 

GTGCUGCTTG TTGAACAACG GTGACCAATT TGGqjGTCTA CGCCGTATCT AG! 




■ I ■ 



■ I 



.1. 



.1 



.1. 



■ I 



.1. 



■ I 



. I . 



190 200 210 220 230 240 

TGGACCGTCG GTAGCGATGA AAGATGCGGA GGAAGTTCCT dT|rcfGj3AGTT GAGCAGACTT 

TGGACCGTCG GTAGCGATGA AAGATGCGGA GGAAGTTCCT CIKG 3AGTT GAGCAGACTT 

TGGACCGTCG GTAGCGATGA AAGATGCGGA GGAAGTTCCT C T TC G 3AGTT GAGCAGACTT 

TGGACCGTCG GTAGCGATGA AAGATGCGGA GGAAGTTCCT C G TC A 3AGTT GAGCAGACTT 

TGGACCGTCG GTAGCGATGA AAGATGCGGA GGAAGTTCCT C G TC A 3AGTT GAGCAGACTT 

TGGACCGTCG GTAGCGATGA AAGATGCGGA GGAAGTTCCT CIKG 3AGTT GAGCAGACTT 

TGGACCGTCG GTAGCGATGA AAGATGCGGA GGAAGTTCCT C G TC A 3AGTT GAGCAGACTT 



typica from dolphin Stenella longirostris, Brasil 
AB479120 A. typica from dolphin Steno bredanensis, USA 
KC928261 A. typica (I) frigate tuna Auxis thazard, Present study 
KC928262 A. typica (II) skipjack tuna Katsuwonus pel amis/ frig ate tuna A. 

thazard, Present study 
EU346093 A. typica from bullet tuna Auxis rochei rochei, western Indonesia 
EU346092 A. typica from bullet tuna A. rochei rochei, western Indonesia 
EU346091 A. typica from bullet tuna A. rochei rochei, western Indonesia 

Fig. 2. Comparison of ITS1 sequences of Anisakis typica larvae in 
the present study. KC928261 A typica (I), KC928262 A typica (II), 
and A typica larva from the bullet tuna from western Indonesia 
(GenBank accession no. EU346093, EU346092, and EU346091 , 
respectively) were compared with adult specimens from dolphins 
from Brazil (AY826724) and USA (AB479120). Square shows dif- 
ferences of nucleotides in the ITS1 regions among A typica. 

on PCR-RFLP and sequencing, Anisakis larvae were identified 
as A. typica. PCR-RFLP was further used to detect proportion of 
Anisakis spp. in fish, and the results showed that from 73 Ani- 
sakis type I specimens collected from K. pelamis (40 specimens) 
and A. thazard (33 specimens), 100% belonged to A. typica. 
From 20 samples of Anisakis type I sequenced, 2 different pat- 
terns, or genotypes, were noted and further identified as A. typ- 
ica. The 2 patterns differed in 4 base pairs in the ITS1 region. 
The first pattern showed 100% similarity with adult A. typica 
(AB479120) reported from dolphins in USA and high similar- 
ity with the Brazilian A. typica (AY826724) from dolphins, 
whereas the second genotype has 4 base pairs different in ITS1 
with A. typica from USA (AB479120), and 100% similarity 
with Indonesian A. typica EU346093 from Auxis rochei rochei. 
The first genotype was only found in A. thazard, whereas the 
second one was found in both K. pelamis and A. thazard. 

Palm et al. [25] proposed that difference in 4 bp of nucleo- 
tide in ITS 1 region may indicate the occunence of A. typica sib- 
ling species. Therefore, the present finding suggests that A. typi- 
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AB551660 J 4. typica- Trichiurus sp 
AB432908 A. typica- Scomber japonicus 
AB479120 A. typica- Steno bredanensis 
AB432909 A. typica —Rastrelliger kanagurta 
EU346092 A. typica- Auxis rochei rochei 
EU346091 A. typica- A. rochei rochei 
EU346093 A. typica-A. rochei rochei 
EU327687 A. typica- Solatia Jluviatilis 

- EU327688A typica- S. Jluviatilis 
EU327686/1. typica- S. Jluviatilis 
EU327689 A. typica- Auxis thazard 
KC928262 A. typica-Katsuwonus pelamis* 
AM706345 A. typica- Selar crumenophthalmus 
KC928261 A. typica-A. thazard* 
AB429051 A. typica 

- FJ234210 A. typica. —Lutjanus lutjanus 
AY826724 A. typica- Stenella longirostris 

AY826725 A. ziphidarum- Ziphius cavirostris 
AY821736 A. simplex C- Mirounga angustirostris 
AY826722 A. simplex C- Pseudorca crassidens 
AY821739 A. simplex C- M. angustirostris 
EU873166 A. pegreJJii-se'Sbre&m 
AB277822 A. simplex- Pleurogrammus azonus 
AF41 1202 A. simplex- Clupae harengus 
AJ937671 A. simplex- C. harengus 
AJ937669 A. simplex- C. harengus 
AJ937670 A. simplex- C. harengus 
-\ EU624343 A.pegreJJii Theragra chalcogramma 
AB277823 A. pegrejjii. S. japonicus 
EU933997 A.pegrejfii. Seriola dumerili 
EW33998 A.pegrejffii 

AM706346 A.pegreJJii Astroconger myriaster 
AY826723 A. simplex. Globicephala melaena 
EU624342 A. simplex s.s. T. chalcogramma 
99 1 AB201 789 A. physteris- pacific cod 

L EU624345 Anisakis paggiae-T. chalcogramma 
r AY826719 A. brevispiculata- Kogia breviceps 
99 ' EU624344 A. brevispiculata- T. chalcogramma 
— EU327690 A. physteris- A. thazard 
AY826721 A. physteris- Physeter catodon 
AB277821 A. physteris-S. japonicus 
L AY603530.4. physteris 

AB211H24-Pseudoterranova decipiens 
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Fig. 3. Phylogenetic tree of Anisakis spp. including the present 
samples (KC928262 A. typica from Katsuwonus pelamis and 
KC928261 A. typica from Auxis thazard) based on ITS-5.8S-ITS2 
gene sequences. Asterisks represent present samples. Maximum 
likelihood tree was constructed using MEGA version 5.1 [33], 
drawn using Nearest-Neighbor-Interchange (NNI) with Kimura 
2-parameter model and 1 ,000 boostrap number with complete 
deletion. Percentages >50% are shown at the internal nodes. 

ca sibling species may occur in K. pelamis and A. thazard. Phy- 
logenetic trees from the mtDNA cox2 region also showed a 
cluster within A. typica, and 2 subgroups were noted. The first 
subgroup consists of all samples from A. thazard and A. typica 
DQ116427, AB517571, and AB517572. The second subgroup 
consists of 4 samples from K. pelamis and 1 sample from A. 
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Table 3. Similarity of nucleotide sequences among A typica including the present Anisakis based on rmtDNA cox2 region 
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Note: 1 -C to 4-C means sample codes of Anisakis typica from Katsuwonus pelamis, and 6-T to 69-T means sample codes of A typica from Auxis 
thazard. 



KC928272 (69-T)* 
DQ116427A fyp/ca-Delphinidac 

- KC928267 (6-T)* 

KC928268 (7-T)* 

— AB517572^(. typica 



- AB517571 A. typica- Scomber 
japonicus 



KC928270 (9-T)* 
- KC928271 (59-T)* 
KC928263 (1-C)* 
KC928264 (2-C)* 
I — KC928265 (3-Q* 
jot— KC928266 (4-Q* 
KC928269 (8-T)* 
DQ1 16430 A. Ziphidanim-ZiphMac 

rDQl 16431 Anisakis sp.-Ziphiidae 
r 




100 L GQl 18171 Anisakis nascettii-Mesoplodon 

grayi -Ziphiidae 

DQ1 16429 A. simplex C-Delphinidae 
DQ1 16428 A. jOegre#;"-Delphinidae 
AJ132189/1. simplex 

DQ1 16426^. simplex s.s- Delphinidae, Phocoenidac, 
Balaenopteridae 

— AF179916 Hysterothylacium reliquens 



Fig. 4. Phylogenetic tree of Anisakis species from the present study 
(KC928263 to KC928272) and other Anisakis spp. based on mtD- 
NA cox2 gene sequences. Asterisks represent present samples. 
Neighbour joining tree was constructed using MEGA version 5.1 
[33], drawn using Maximum composite likelihood Model and 1 ,000 
boostrap number with complete deletion. Percentages >50% are 
shown at the internal nodes. Sample codes were presented in Ta- 
ble 3. 

thazard. The first subgroup of the present samples showed 96- 
100% similarity with the known A. typica, while the second 
subgroup had 93-95% similarity. The second subgroup forms 
another cluster which separate them from the known A. typica 
and this may also indicate the presence of A. typica sibling spe- 



cies. A previous study in Papua New Guinean waters also 
showed a similar pattern with the present study in which ge- 
netic divergence occuned within A. typica clade [34]. Though a 
study in other nematode taxa showed that sequence differenc- 
es of about 10-20% were interspecific, and differences of about 
7% were regarded as conspecific [35], the present study, as 
proposed by Palm et al. [25], indicates the presence of A. typica 
sibling species in K. pelamis and A. thazard, whereas A. typica 
was recorded only from A. thazard. 

Molecular differentiation of Anisakis spp. using PCR-RFLP 
has been successfully used [16-20,27,35-38]. In the present 
study, digestion of 73 samples with restriction enzyme Taq I, 
Hinfl, and Cfo I indicated that the samples were A. typica, and 
suggested a predominance of A. typica in the eastern parts of 
Indonesia. This high abundance of A. typica in the present study, 
and the report of A. typica in the previous study in Balinese 
and Javanese waters [25], as well as a recent report from Papua 
New Guinea [34] support the previous findings that this spe- 
cies was more abundant in temperate and tropical waters [39]. 
Mattiucci et al. [39] recorded larvae of A. typica from A. thazard 
and Thunnus thymus from Brazilian Atlantic Ocean, Scomber 
japonicus and Trachurus picturatus from Madeira Atlantic Ocean, 
E. affinis, S. commerson, Sarda orientalis, and Coryphaena hippu- 
rus from Somali Coast, and Merluccius merluccius from the 
Mediteranian Sea off Cyprus and off Crete. In Indonesia, more 
than 34 species of marine fish have been reported to harbor 
Anisakis spp. [25]. Recent reports of Anisakis infection added 3 
more new fish genera to harbor Anisakis spp. Two genera were 
reported from the Southern coast of Kulon Progo; Parupeneus 
sp. (Mullidae) and Terapon jarbua (Terapontidae) [26], and 1 
genus in the present study, C. cyanostigma (Senanidae), is a new 
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record of Anisakis infection. At present, based on molecular 
studies, 9 species of Anisakis are known, namely A. simplex s. s., 
A. pegreffii, A. simplex C, A. typica, A. ziphidarum, A. nascettii, A. 
physeteris, A. brevispiculata, and A. paggiae [21,22]. 

A. typica has been reported from numerous marine fish 
worldwide. The parasite was reported from marine fish in Ko- 
rea, Japan, China, Portugal, Taiwan, Brazil, Western Indonesia, 
Morocco, Papua New Guinea, Adriatic Sea of Croatia, Mauri- 
tania, and some countries at Mediterranean Sea [20,25,34-37, 
40-47]. The existence of A. typica from the Portuguese coast may 
have extended the distribution of this parasite to cold water. 
However, the infection level of the parasite was very low. There- 
fore, it might be possible that the fish may accidently infected 
through the food chain originating from warm waters. Marques 
et al. [45] stated that the Portuguese coast is a transition be- 
tween North-Eastern Atlantic warmer temperate and cold tem- 
perate regions so that it might provide an area of species over- 
lap and hence could promote accidental infection. 

High prevalence of Anisakis was found in migratory fish, K. 
pelamis and A. thazard. High prevalence of Anisakis infection 
was also reported from some marine fish in southern coast of 
Kulon Progo, Yogyakarta [26]. The prevalence of Anisakis sp. 
infection was generally higher in bigger fish than in smaller 
ones [19,25,48]. The same result was found in this study that 
small A. thazard were not infected with Anisakis, while bigger 
ones were infected with the prevalence of 47% and the mean 
intensity of 5.6. This result might be caused by accumulation 
of the parasites in the big fish due to a long period of infection. 
Previous reports on Anisakis infection in Indonesian waters 
showed high prevalence of infection with the parasite in some 
species of marine fish. Hadidjaja et al. [24] reported that the 
prevalence of Anisakis type I larvae in Rastrelliger kanagurta, De- 
capterus russeli, and Sardinella sirm was 49.3%, 50.3%, and 
40.9%, respectively, whereas in the present study no Anisakis 
infection was found in D. russeli and only 5% infection in ft. 
kanagurta. The difference in the prevalence of infection was 
also noticed at different locations by Palm et al. [25], and they 
suggested that the high prevalence of Anisakis infection at 
Northern Balinese coast was due to the high abundance of 
dolphins, as the final host for Anisakis, in that area. A previous 
study on the ecology of Pseudotenanova decipiens in Antarctic 
waters showed that a high prevalence of infection was in ac- 
cordance with a high abundance of final hosts as well as inter- 
mediate hosts in the area [49]. 

Anisakiasis has been reported from several countries such as 



Japan, Korea, and some countries in Europe. Anisakis may in- 
fect humans and causes anisakiasis after consuming raw in- 
fected fish or other marine organisms that function as inter- 
mediate hosts. The first report of anisakiasis was from a patient 
in the Netherlands who had gastrointestinal problems due to 
A. simplex infection. Most cases of anisakiasis in Europe and 
Japan have been reported to be caused by Anisakis type I, par- 
ticularly A. simplex [4]. However, anisakiasis due to A. pegreffii 
infection has also been reported from humans in Italy [50,51]. 
A. simplex might penetrate and migrate to fish muscle [52], 
which may explain the higher cases of anisakiasis due to A. 
simplex infection. Anisakiasis due to A. typica has not been re- 
ported. Umehara et al. [20] reported that A. typica so far has 
only received limited attention and is not widely recognized, 
thus its zoonotic impact has not been well documented. How- 
ever, Palm et al. [25] reported that A. typica was not only found 
on the surface of gastrointestinal ttact but it might also pene- 
trate muscle of fish. In the present study, though not common, 
an Anisakis larva was observed to migrate into the musculature, 
indicating that the parasite has the potential to infect humans 
through consumption of uncooked food. In Indonesia, reports 
about anisakiasis in humans were suggested by Uga et al. [9] 
using a seroepidemiological approach of inhabitants in East 
Java and revealed that about 11% of samples who visited hos- 
pital showed positive results for Anisakis antibodies. The spe- 
cies of Anisakis spp. was not determined, but it might be possi- 
ble that the parasite was A. typica since this species is widely 
distributed in tropical waters, compared to the well known 
causative agent of anisakiasis by A. simplex and A. pegreffii which 
have limited geographical distribution in cold waters. Howev- 
er, it might also be possible that Anisakis spp. could be from 
imported raw materials. Yoshinaga et al. [53] reported the 
presence of A. pegreffii in amberjack Seriola dumerili imported 
from China to Japan as mariculture seedlings. 

In conclusion, based on PCR-RFLP and sequencing, all the 
Anisakis examined were A. typica, indicating the predominance 
of this species in the Southern Makassar Strait, Indonesia. Se- 
quencing and phylogenetic tree analyses of Anisakis type I in 
ITS1-5.8S-ITS2 and mtDNA cox2 regions showed that the pres- 
ent samples were in the same cluster as A. typica published in 
GenBank. However, differences of 4 bp in nucleotides in ITS1 
region and broad divergence consisting of 2 subgroups in the 
mtDNA coxl of Anisakis from K. pelamis and A. thazard indicat- 
ed the existence of A. typica sibling species in that area. 
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